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•  The  West  Desert  Test  Center  (WDTC)  at  DPG 
uses  an  elastic  backscatter  LIDAR  system  to 
augment  the  tracking  and  characterization  of  bio¬ 
aerosol  clouds. 

•  In  support  of  DPG’ s  test  mission,  hardware  and 
software  development  for  WDTC  LIDAR 
capability  is  continuing. 

•  This  talk  is  a  progress  report  on  our  M&S  project 
to  estimate  bio-aerosol  concentration  from  LIDAR 
elastic  backscatter  data 
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A  convenient  form  of  the  Lidar  equation  pre¬ 
dicting  the  received  power  (in  Watts),  from  an 
increment  of  atmosphere  A R  at  range  R,  is 


P(A,  R)  =  P0 


Here  A  is  wave  length,  Aq  is  area  of  the  receiver 
telescope  objective  lens,  fa  is  the  backscatter 
coefficient  and  ne  is  the  extinction  coefficient. 
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Backscatter  Coefficient 


The  backscatter  coefficient  (5-n  (m  1  sr  x)  is 

r  oo 

/?tt(  A,  r)  =  /  iV(a,  r)  cr^(A;  a)  da  . 

^  o 

Distribution  iV(a,  r)  is  normalized  such  that 

roc 

p(r)  =  /  7V(a,  r)  da  . 

JO 

Here  p(r)  is  number  of  particles  per  unit  vol¬ 
ume.  The  Mie  backscattering  cross-section  is 


<t6(A,  a) 


1  .  1 


OO 


fc2  1  2 


-  XI  (2m+l)(-l)m(am-bm) 
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m=  1 
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The  extinction  coefficient  /ce(A,r)  (m  *)  is 


roc 

/ce(A,r)  =  p(r)  /  Ni(a)ae(\a)da 

Jo 


where  the  extinction  cross  section  is 


2tt  00 

(Te(\o)  =  (2m  +  1)  {dm  +  bm)  ]  • 

l\j  _ -i 

m=  1 
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Inversion  of  LIDAR  Data 


To  estimate  aerosol  concentration,  transform 
the  measured  LIDAR  data  to  obtain 


S(R,  A)  =  log(/?2JP(A,  R)/(R%R(\,  R0 ))  . 
This  definition  then  yields  the  explicit  relation 


-j—S(R,  A)  -  — -  2  p(i?)  <  cre( A)  >  . 
art  p{lt)  art 

This  equation  is  solved  for  the  particle  density 
function  p( r).  Necessary  information  includes 
the  averaged  extinction  coefficient  <  cre(A)  >, 
LIDAR  data  S(R,  A)  and  a  boundary  point  pf. 


The  Nation’s  Chemical  &  Biological  Defense  Proving  Ground 


Formal  Analytic  Solution  to 
Aerosol  Density 


tRt  dE1  i 

p(R)  =  [(l/p/4-2  <  oe(\)  >  Jr  — ) r(R,  A)]"1 , 

where 


T 


The  Nation’s  Chemical  &  Biological  Defense  Proving  Ground 


Aerosol  Scattering 


Mie  Efficiencies  for  n  =  1.33 


08 


Lidar  signals  are  proportional  to  Qb  and  exponentially 
attenuated  proportional  to  Q 
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Aerosol  Particle  Size 
Distribution 


a,  particle  radius  (n  m),  N  =100  cm  3  02-^-2002  aqh 
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Extinction  Cross  Section  versus 
Aerosol  Mode  Radius  a. 


m 
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Particle  Density  Function 


Gaussian  Cloud  Density  p(  r)  PPL 
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Synthetic  Density  Data 
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LIDAR  Data 


Free  Lidar  Data  S 
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First  Stage  Inversion 
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First  stage  inversion  results  using  two-jrequency  method. 


Second  Stage  Inversion 
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Second  Stage  inversion  with  S/N  =  16  iS  noise  using  eguations  for  (jnm  and  point  detector 
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Linear  Scale  Inversion  Results 


Particle  Density  (particles/liter) 
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Estimated  Particle  Density,  S/N  =  10  dB 


Range  R  (m)  «j>  (Deg)  30^Apr^2003  aqh 


Linear  scale  particle  density  inversion  in  S/N  =  10  dB  noise.  Upper  sub-plot  is  exact  density. 


Two  Frequency  Method 


Reference  Arc  Estimates 
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Figure  11:  Inversion  error  nnnlysis  ns  n  function  of  nzimuth  nngle  frnction  with 

five  vnlues  of  ns  shown  in  the  legend .  Upper  left-hnnd  pnnel  is  the  relntive  percent 
error  in  the  extinction  coefficient  estimnte  o^1)  =  2  <  rre( A  i  )  >.  The  upper  right  pnnel 
displnys  the  nnnlogous  percent  relntive  errors  in  nvernged  honndnry  vnlues  q f .  Note 
the  condition  number  is  independent  of  /A . 
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Discrete  Forward  Model 


where 


Density  Calibration 
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Single-Ray,  Two-Frequency 
Algorithm 
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Single-Ray,  Three-Frequency  Algorithm 
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Conclusion  &  Future  Work 


•  2  and  3  frequency  elastic  scattering  algorithms 
for  particle  density  have  been  developed. 

•  At  least  2  frequencies  are  necessary  to 
determine  reference  arc  values  pf. 

•  Ongoing  work  includes  background  removal, 
first  stage  signal-to-noise  enhancement,  and 
system  conditioning. 
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